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j^ , Abstract I review the basic physics of ultracold dilute trapped atomic gases, with 

bX)' emphasis on Bose-Einstein condensation and quantized vortices. The hydrody- 

namic form of the Gross-Pitaevskii equation (a nonlinear Schrodinger equation) 
P5 ■ illuminates the role of the density and the quantum-mechanical phase. One unique 

C^ ' feature of these experimental systems is the opportunity to study the dynamics 

^ . of vortices in real time, in contrast to typical experiments on superfluid ^He. I 

discuss three specific examples (precession of single vortices, motion of vortex 
"ti [ dipoles, and Tkachenko oscillations of a vortex array). Other unusual features in- 

clude the study of quantum turbulence and the behavior for rapid rotation, when 
the vortices form dense regular arrays. Ultimately, the system is predicted to make 
1^ . a quantum phase transition to various highly correlated many-body states (anal- 

P3 ' ogous to bosonic quantum Hall states) that are not superfluid and do not have 

O , condensate wave functions. At present, this transition remains elusive. Conceiv- 

ably, laser-induced synthetic vector potentials can serve to reach this intriguing 
phase transition. 

PACS numbers: 03.75.Hh, 05.30.Jp, 67.40.Db 



1 Introduction 

Consider a uniform gas of particles with mass M and number density n. The inter- 
particle spacing is ^ n^''-^, and there are two different approaches to discuss the 



^^ ' onset of quantum degeneracy (for general background, se c^'^i^'^i^ ). 

Start with an atomic -physics perspective: the mean thermal energy p^ /2M : 



kgT yields a mean thermal momentum p ^ ^/MkgT, and the familiar de Broglie 
relation X ^ h/ p gives the mean thermal wavelength A ^ h / \/ MkgT . Compare 
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A with the interparticle spacing n ''^. In the classical limit (high temperature, 
short wavelength), we have A ^ rT^'^ . Hence quantum diffraction is negligible, 
similar to ray optics for light when diffraction is unimportant. It is convenient 
to define the dimensionless parameter nX^ , known as the "phase-space density;" 
this parameter is small in the classical limit since A — ;> when h^O ov T ^ °o. 
As the temperature falls at fixed n, the thermal wavelength grows, and quantum 
degeneracy appears when nk^ is of order 1. Forbosonic atoms, this criterion yields 
the temperature Tc for the onset of Bose-Einstein condensation. For T < Tc, a 
macroscopic number of particles No{T) occupies the lowest single-particle state, 
and the fraction of particles in this lowest state increases as T decreases. In an ideal 
gas, all particles occupy the single-particle ground state at T = K. Dilute trapped 
quantum gases have low densities with n ^ 10^^ cm~^, roughly 10"^ smaller than 
room-temperature air. This density gives a low transition temperature Tc ^ I0~^ K. 

Alternatively, take a condensed-matter view. Each particle occupies a "box" 
of dimension n~^'^, which gives a zero-point confinement energy ^ h rP''^ /M. In 
the classical limit, the thermal energy ksT is much larger than the confinement 
energy, but as T falls, the system eventually reaches the transition temperature 
ksTc ^ trn^l^ /M for the onset of quantum degeneracy. For a Bose system, this 
analysis gives the same criterion for Tc as found from the phase-space density. For 
a Fermi system, it gives the usual Fermi temperature Tf. Electrons in metals have 
a large Tf ^ 10'* K, whereas liquid ^He has Tf ^ \}L because the number density 
remains similar but the mass is larger by roughly 10'*. Typical dilute trapped Fermi 
gases have low number density and large mass, which leads to Tf ^ 10~^ K, like 
a dilute trapped Bose gas. 

These ideas apply directly to an ideal Bose gas in a spherical harmonic trap, 
with trap potential V^{r) = ^May^r^. The familiar ground-state wave function is 
a Gaussian V'b('') °^ sxp{—r^/2d^), with the characteristic size given by the os- 
cillator length d = y^h/{M(o). For dilute atomic gases, d is typically a few /im. 
In a harmonic trap with A'^ particles, the onset of BEC occurs at kgTc ^ hcoN^'^; 
for the typical value N ^ 10^, this yields Tc ^ I /iK. Above Tc there is only a 
wide thermal cloud, but below Tc a narrow condensate of width d starts to ap- 
pear, rising from the much wider thermal cloud. As T — ;> K, the thermal cloud 
disappears, leaving only the narrow condensate. The presence of a Bose-Einstein 
condensate at T = K means that all the A'^ condensed particles form a coher- 
ent quantum state described by a macroscopic wave function f (r) = \/N\j/Q{r). 
The normalization is J dV |f p = N for T <C Tc. In this ideal gas, the condensate 
density (which is the same as the total density at low temperature) is nonuniform 
with n{r) = |f (r)p. In a spherical harmonic trap, the condensate density has the 
Gaussian form n{r) <== NQX\){—r^/d^). 

What is the effect of interparticle interactions? The basic idea was due to Bo- 
goliubov: for weak interparticle potentials and T ^ Tc, nearly all the particles 
remain in the condensate. For cold dilute gases, the typical interparticle spac- 
ing is n~^'^ ^ a few 100 nm. In contrast, the interactions are short range, char- 
acterized by the i-wave scattering length a ^ a few nm. Thus the dimension- 
less parameter no? is small (of order 10~^ for a typical cold dilute gas). Each 
particle in the condensate experiences a Hartree mean-field potential Vh(''") = 
{4nah^/M)n{r) = {4nah^/M)\W\^ from the contact interaction with all the re- 



maining condensed particles. The self-consistent condensate wave function obeys 
a nonlinear Schrodinger equation, usually called the Gross-Pitaevskii (GP) equa- 
tion 

where fi is the chemical potential. 

In the presence of a trap, this GP equation involves a new dimensionless "in- 
teraction" parameter Na/d that combines the effects of the trap and the interac- 
tions. Note that the typical ratio a/d is small, of order 10"^, but with A' w 10 , 
this GP interaction parameter is large. Hence the repulsive interactions expand the 
condensate to a radius R that significantly exceeds the ideal radius d (a typical 
value is R/d ^ 10). Consequently, the radial gradient of f becomes small when 
Na/d ^ 1 . In this limit, one can neglect the kinetic energy term in the GP equation 
dlj,— leading to the much simpler algebraic relation 

l^|f(r)|2=^-y,(r), (2) 

which is called the Thomas-Fermi (TF) approximation. For a spherical harmonic 
trap, this result yields the simple particle density (an inverted parabola) 

n(r)=«(0)(^l--^), (3) 

where n(0) = Mii/{Anah^) is the central density and R^ = 2pL/{Map-) is the 
squared condensate radius {R is the classical turning point for a particle with en- 
ergy ii). 

Section 2 reviews the time-dependent Gross-Pitaevskii equation, including the 
structure of a single vortex in an unbounded condensate. The more complicated 
situation of a single vortex in a large trapped condensate is analyzed in Sec. 3, 
and Sec. 4 summarizes some of the experiments on creation and detection of vor- 
tices, including large arrays. Turbulent vortex systems in trapped condensates are 
considered in Sec. 5, where at present only a few experiments exist. Sections 6 
and 7 focus on vortex arrays in the mean-field regime, with two distinct cases: the 
Thomas-Fermi regime when the vortex cores are well separated and the density 
variation is negligible, and the lowest Landau level regime, when the vortex cores 
overlap and the spatial variation of the density becomes important. In principle, 
the system should make a quantum-phase transition to a highly correlated non- 
superfluid state for sufficiently fast rotation (Sec. 8), but experiments have not yet 
achieved this limit. 



2 Time-dependent Gross-Pitaevskii equation 

This central nonlinear equation was proposed independently by Gross and Pitaevskii 
in 1961.^'^ It provides an accurate description of low-temperature trapped Bose- 
Einstein condensates, ^-^ in particular 

1 . properties of the ground state. 



2. free expansion of the condensate after the confining trap is turned off, 

3. collapse for attractive interactions, 

4. frequency of low-lying collective modes (at the 1% level of accuracy). 

It comes in two different but equivalent versions. 



2.1 Nonlinear Schrodinger equation 

This view emphasizes the quantum aspects of the problem 






It can be considered to arise from an energy functional 



E[w]=ldV[^^+V,\W\^ + ^\Wr]. (5) 



The equilibrium state minimizes this energy functional with fixed normalization 
J dV |f 1^ = A'^, yielding Eq. ([U as the resulting Euler-Lagrange equation. In this 
approach, the chemical /i in Eq. ([TJ serves as a Lagrange multiplier that enforces 
the constraint of fixed A'^. 



2.2 Hydrodynamic description 

Write the condensate wave function as f = exp(i5')|f |, which yields hydrody- 
namic variables: 

1. particle density n = \W\^, 

2. velocity v = hVS/M, which is irrotational except for singularities, 

3. particle current density j = nv. 

The circulation is defined as K = §^dl ■ v on a closed path C As in superfluid 
^He, the single-valued condensate wave function implies quantized circulation 
with k: = integer x 2nh/M (note a similar result holds for type-II superconductors, 
where the electronic charge plays a crucial role). 

Substitute the form f = exp(j\S') |f | into the time-dependent GP equation Q. 
The imaginary part gives the usual conservation of particles dn/dt + V • {nv) — 
0. The real part gives a generalized Bernoulli equation that incorporates all the 
physics of compressible irrotational isentropic hydrodynamics, including quan- 
tized vortices and their dynamics. The one new feature is that the trap potential Vtr 
makes an additional contribution to the vortex motion. 



2.3 Straight singly quantized vortex in bulk fluid 

Gross ^ and Pitaevskii^ independently studied the structure of a singly quantized 
vortex in an unbounded condensate with bulk density n. This weakly interacting 
dilute gas provided a toy model for a vortex in dense superfluid ^He. Assume a 
condensate wave function f (r) = ^/nexp{i^)f{r), where r is in the xy plane with 
^ the polar angle. This choice yields a flow with circular streamlines 

h - 

Note that v diverges as r — !■ 0. The circulation is (c = 2nh/M with singular vor- 
ticity at the origin V x t; = JciS'^^ (r). The chemical potential is ^u = AKatP'n/M 
which is the Hartree energy for a uniform fluid. The speed of sound is simply 
i = Y^/j/M' 3nd real .s requires a repulsive interaction with a > 0. The centrifugal 
barrier forces the radial function /(r) to vanish at the origin with a core radius 
w the healing length ^ fixed by the balance between the kinetic energy and the 
interaction energy jj. 

The vortex core is large compared to the interparticle spacing for a dilute gas with 
na^ <^ 1 . The circulating flow (|6]l around the vortex becomes supersonic near the 
core; from this perspective, the vortex core arises from acoustic cavitation. 

Note that Gross and Pitaevskii sought to model a vortex in superfluid "^He, 
which has only a single superfluid component. For dilute trapped alkalai-metal 
gases like ^Li, ^^Na and **^Rb, however, it is easy to make two-component mix- 
tures with various hyperfine states arising from the spin of the single unpaired 
valence electron and the nuclear spin. In this case, the vortices can have more 
complicated internal structures, along with arrays containing vortices in each com- 
ponent. Section IV.B.5 of Ref.^ briefly summarizes the current situation, in partic- 
ular the dynamical experiments of Schweikhard et al. ^ For simplicity, this article 
will focus on single-component vortices, although the first experimental creation 
of a vortex in a trapped condensate in fact relied on the coupling between two 
hyperfine states in **^Rb (see Fig. 2 below). 



3 Single Vortex in a Large Trapped Condensate 

Assume an axisymmetric harmonic trap potential with 

Vtr(r) = V^{r,z) = \M (ft)2 r^ + (o^) . (8) 

For (0; ^ ft)^, the condensate has a flattened disk shape because of the strong axial 
confinement, whereas for ft), <C (0^, the condensate has an elongated cigar shape 
because of the strong radial confinement. 

The usual experimental procedure is to turn off the trap potential and then take 
an image of the expanded condensate at a later time. For definiteness, consider a 
cigar shape with ft)^ ^ ft),. Before the expansion, the tight radial confinement has 
a large radial potential energy and a small axial potential energy. After the trap is 



turned off, the condensate expands rapidly in the radial direction and soon acquires 
a flattened disk shape. Similarly, an initial disk-shaped condensate expands axially 
to become elongated. This behavior is known as the "reversal of aspect ratio;" it 
provided crucial evidence for the existence of a BEC in early experiments (see^). 
Assume that such a trap rotates around z with angular speed £1. In the rotating 
frame, the original Hamiltonian changes ^'^ from H to H' = H i2L,. For equi- 
librium in this rotating frame, the condensate wave function W minimizes the GP 
energy functional in the rotating frame 



E'[W,Q]=E[W]- fdVW*{f2-rxp)W. 



(9) 



For simplicity, assume a rotating disk-shaped condensate of radius R±. Use E' = 
E — D.L, to determine the energy Eg of a rotating condensate with no vortex and 
the energy E[ of a rotating condensate with a straight off-center vortex at a distance 
tq from the center. Define the creation energy for the vortex as Zii?' = iij — ZJq. For 
a typical dilute gas, the core radius ^ is a few x 10~^ m. Detailed analysis ^' yields 
the following picture (Fig. 1) of the vortex energy AE' as a function of the radial 
displacement tq for various fixed values of Q. . 



A£'(roJl) 




Fig. 1 Energy of a single off-center vortex in a disk-shaped condensate for various fixed 12: (a) 
CI =0,{b) Q. = Q.i„ for onset of metastability, (c) .Q = 12^ for onset of tiiermodynamic stability— 
(reprinted with permission of the author and the American Physical Society). 



Curve (a) is for 12 = 0, when the energy AE' decreases monotonically with in- 
creasing ro, and the trap center is a local maximum of this curve. In the absence of 
dissipation, a fixed energy means a fixed radial position tq. Hence the only allowed 
motion for such a vortex is uniform circular precession. The rate of precession 
is proportional to the slope of the energy curve at tq with ^ oc —dAE'/drQ. In the 
presence of weak dissipation, the vortex moves down the energy curve to reduce 
its energy, so that the vortex slowly spirals out of the nonrotating condensate. 



With increasing external Q, the negative central curvature of the energy curve 
in Fig. 1 decreases, and curve (b) is the special value 

n,„ = l_^i^(^] (10) 

when the central curvature vanishes. This result means that a central vortex first 
becomes metastable at I2„,, because the trap center is a local minimum of the en- 
ergy for Q > Q,„. For weak dissipation, a vortex would now spiral inward for 
small lateral displacements to reduce its energy. Note that such a vortex is not 
globally stable since AE' remains positive at the origin. Curve (c) occurs for 
Qc = f An; when the vortex first becomes truly stable. 

How fast does an off-center vortex precess? One convenient approach relies 
on the Lagrangian variational approach, ^^i^^ where the Lagrangian functional ^ 
has the detailed form 

if, f / dW dW* \ 

This functional is stationary for small variations of W and f *, with the familiar 
time-dependent Gross-Pitaevskii equation as the Euler-Lagrange equation. If the 
trial wave function contains one or more parameters, the resulting ^ contains the 
time derivative of the parameters because of the explicit appearance of dW/dt and 
dW* /dt. The resulting Lagrangian serves to study the dynamical evolution of the 
parameters. 

In a disk-shaped TF condensate with radius R^, the vortex position ro(f ) pro- 
vides a simple example of such a time-dependent parameter. For a nonrotating 
condensate, the precession rate becomes 



where i2,„ = j(Ji/MR^)ln{R^/^) is the frequency for the onset of metastability 
given in Eq. (llOb . It is notable that the precession is positive, in the same sense as 
the fluid flow around the core. Here, the factor 1 — r'^/R\ comes from the parabolic 
radial TF density profile (ultimately this dependence reflects the quadratic har- 
monic trap). As discussed below, experiments on trapped BECs confirm this re- 
sult in considerable detail. It is instructive to compare this result with that for a 
similar vortex in incompressible fluid bounded by a rigid cylinder of radius /?^ (as 
a model for superfluid ^He). Classical hydrodynamics yields a formally similar 
result 

but here the denominator arises from the image vortex located at a distance R\/rQ, 
instead of the nonuniform density. Also, there is no "large" logarithmic factor 

HRJ^). 



4 Experimental creation and detection of vortices in BECs 

Most such experiments study equilibrium vortex arrays. The first vortex was made 
at JILA (Boulder, CO) in 1999. ^^ They used a nearly spherical **^Rb condensate 
containing two different hyperfine states. A coherent laser coupling controlled the 
interconversion between the two species, and a stirring perturbation could spin up 
the condensate. When the laser coupling was turned off, they obtained one com- 
ponent with a singly quantized vortex that circulated around a nonrotating core 
of the other component. Selective laser tuning provided nondestructive images of 
either component. These images (Fig. 2) allowed a study of the precession of such 
a two-component vortex around the trap center. '^ The core can be as large as 5- 
10 /J m, depending on the fraction of the nonrotating component, which is readily 
imaged with visible light. 
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Fig. 2 Precessing two-component vortex (a) direct images, at 50 ms intervals (b) smoothed 
images, (c) angular position, and (d) gradual shrinkage of core due to spin-flip transitionsi^ 
(reprinted with permission of the authors and the American Physical Society). 



The JILA group could also remove the nonrotating core component with an 
intense laser pulse, leaving a single component vortex with an empty core, which 
marked the initial position. They waited a variable time and then turned off the 
trap, allowing a visualization of the final position of the vortex. The observed 
precession rate agreed with the theoretical analysis at the ±10% level. The exper- 
iments saw no outward radial motion for '^ 1 s, implying that dissipation is small 
on this time scale. 

The Ecole Normale Superieure (ENS) group in Paris studied vortex creation in 
a very elongated rotating cigar-shaped condensate with one component. They used 
an off-center toggled rotating laser beam to deform the transverse trap potential 



□□ 



T # y 



Fig. 3 Images of small vortex clusters in a rotating BEGiS (reprinted with permission of the 
authors and the American Physical Society). 



and stir the condensate at an applied frequency Q./2k < 200 Hz. They observed 
small arrays of up to 11 vortices arranged in two concentric circles (Fig. 3). They 
needed to expand the condensate to obtain these pictures. ^^ Note that these images 
are like patterns predicted and seen in superfluid "*He (Fig. 4).— 




Fig. 4 Images of small vortex clusters in rotating superfluid *He " (reprinted with permission 
of the authors and the American Physical Society). 



Soon afterward, the MIT group prepared considerably larger rotating conden- 
sates in a less elongated trap. '^ They observed large triangular arrays with up to 
130 vortices, like the Abrikosov vortices of quantized flux lines in type-II super- 
conductors. As an alternative approach, Cornell's group at JILA started from a 
rapidly rotating normal cloud and then cooled into the superfluid state that con- 
tained a vortex array to accommodate the large angular momentum (Fig. 5).— 122 

Anderson's group in Arizona has created vortex dipoles (± vortex pairs) in a 
disk-shaped TF condensate.^' They used an intense blue-detuned laser beam as 
an obstacle in the condensate, moving the condensate so that it sweeps smoothly 
past the obstacle at differing rates. Above a critical rate, they created a vortex 
dipole at a reproducible position, waited a variable time, and then turned off the 
trap. Figure 6 shows expanded pictures of the dynamical motion of the vortex 
dipole at intervals of 200 ms. The first row is experimental data, and the second 
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Fig. 5 Large triangular vortex array in a rotating BEC^° (reprinted with permission of the au- 
thors and the American Physical Society). 
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Fig. 6 Creation of vortex dipole by a blue-detuned laser beam in a disk-shaped TF condensate, 
(a) Experimental data, (b) numerical simulation, (c) comparison of observed and theoretical 
trajectory— (reprinted with permission of the authors and the American Physical Society). 



row is a theoretical simulation. The figure (c) on the right compares the measured 
trajectory (dots) with the theoretical trajectory (continuous line) 



5 'Rirbulent vortex systems 

Kobayashi and Tsubota^^ have carried out numerical simulations of the time- 
dependent GP equation with rotations first about the z axis and then about the 
X axis, leading to an effectively time-dependent rotation axis f2{t) because of the 
combined rotations. For a slightly asymmetric triaxial condensate, they find that 
the condensate surface eventually becomes irregular, and a turbulent vortex tangle 
then develops for sufficiently long times. 

Bagnato's group in Sao Carlos, Brazil use a related scheme of oscillations 
about two axes to generate what appears to be a turbulent vortex tangle in a cigar- 
shaped condensate. ^^ When the trap is turned off, the "turbulent" condensate ex- 
pands with an approximately self-similar profile, in contrast to the usual reversal 
of aspect ratio for a nonrotating condensate (Fig. 7). For a uniform vortex array, 
the vorticity induces an additional expansion in the two perpendicular directions. 
Thus a combination of random turbulent vorticity in all three directions may ex- 
plain this unexpected self-similar behavior. 



6 Vortex arrays in mean-field Tliomas-Fermi regime 

As expected, the mean vortex density n^, in rotating condensates obeys the Feyn- 
man relation familiar from superfluid "^He, with n^ = 2X2 /k" = MQ. /{nh). Corre- 
spondingly, the area per vortex is l/n,, = nh/{MQ) = nr, which defines the ra- 
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Fig. 7 Comparison of expansion of nonrotating condensate (left side) and "turbulent" conden- 
sate (right side) ^^ (reprinted with permission of the authors and the American Physical Society). 



dius I — y^h/{MQ) of an equivalent circular cell. Note that the intervortex spacing 
'^ 21 decreases like 1/^/Q. 

With increasing Q, the mean vortex density grows linearly following the Feyn- 
man relation. In addition, centrifugal forces expand the condensate radially, so 
that the area TtR^ also increases. Hence the number of vortices A^,, = M£2R^/h 
increases faster than linearly with Q . Conservation of particles implies that the 
condensate also shrinks axially. The Thomas-Fermi approximation assumes that 
the interaction energy (^^If |^) and the trap energy (Vtrlf P) are both large rel- 
ative to the gradient energy for density variations (/i^/2M)((V|f |)^). This TF 
approximation holds for well-separated vortices with I ^ ^, but it breaks down 
when the vortex lattice becomes "dense" and the cores start to overlap. 

For a quantitative description, note that the kinetic energy involves — /Vf w 
MvW/h from the gradient of the phase, since the density variation is here negligi- 
ble. Hence the TF energy functional in the rotating frame in Eq. (|9]l becomes 

E'[W]= fdV[{\Mv^ + Vt,-Mf2-rxv)\W\^ + ^g\W\'^], (14) 

where g = 4nah /M is the coupling constant and v is the velocity generated by 
all the vortices. In the present limit of many vortices, the Feynman relation im- 
plies that this velocity is just the solid-body result v^\, = fi xr. For Q. along z, 
substitution into Eq. (fT4ll yields 

E'[W]= j dV[\M{(ol-n^)\W\'^ + \M(ol\W\^ + \g\W\'^], (15) 

which now looks exactly like the TF energy for a nonrotating condensate, but with 
a reduced squared radial trap frequency (o\ — > (o\ Q.^. 

Hence the TF condensate density now depends explicitly on Q.: |f(r,z)p = 

n(0) (l - r^/R\-z^/R^), with 



,2 2;U 2 2jU 

~ — ; — ^ ^ and R. = ^ 

M{(ol - 122) z Mcol 



^i = TTTT^I^^T^ ™d Rz = T7^- (16) 



It is clear that Q cannot exceed (O^, since otherwise the radial confinement would 
disappear. In addition, the central density «(0) and the chemical potential jJL — 
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gn{0) both decrease with increasing Q because of the reduced radial confinement. 
The formulas ( fT6l) for the condensate radii show how the aspect ratio changes with 

n 



Rim 



(0, 



122 



(0, 



(17) 



This last effect provides an important diagnostic tool to determine the actual an- 
gular velocity Q (Fig. 8).^"* The measured aspect ratio indicates that Q/o)±^ can 
become as large as w 0.993. 




Fig. 8 Increasing angular velocity dramatically flattens an initially cigar-shaped condensate— 
(reprinted with permission of the authors and the American Physical Society). 



In 1966, Tkachenko^ studied the collective modes of an infinite triangular 
vortex lattice for motion perpendicular to the vortex axes. He predicted a special 
mode involving long-wavelength transverse shearing motion of the straight vortex 
lines. This behavior arises from the discrete quantized vorticity in each vortex. It 
disappears for "inertial waves" in a rotating classical fluid with uniform vorticity 
(namely, k: — ;■ such that «,,k: = 2Q). 

Cornell's group at JILA studied the dynamical motion for such Tkachenko 
waves in trapped BECs. They formed a uniform vortex array and then applied 
a weak perturbation,^^ setting up these Tkachenko waves. Figure 9 shows the 
deformed vortex lattice at ^ and | of the oscillation period, with clear evidence of 
the phase reversal. The observed motion has the correct quantitative form, but the 
measured period differs from the predictions. Sonin's review article^^ discusses 
earlier studies of Tkachenko waves in superfluid ^He (see Sec. VI. E, particularly 
Fig. 3). 




Fig. 9 Deformed vortex lattice at 1/4 and 3/4 period, showing the transverse Tkachenko wave— 
(reprinted with permission of the authors and the American Physical Society). 
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7 Vortex arrays in mean-field lowest Landau level regime 

When the vortex cores begin to overlap, it becomes necessary to include the ki- 
netic energy arising from the density variation near each vortex core. Evidently, 
the Thomas-Fermi approximation breaks down, for it ignores such rapid density 
variations. Hence it is essential to return to the full GP energy functional [Eq. (|9]l] 
for £'[f ,12] in the rotating frame. In this rapid-rotation limit {Q < ft)^), Ho^- 
pointed out that it is possible to incorporate the full kinetic energy exactly. The 
centrifugal forces expand the condensate, which becomes disk-shaped. For sim- 
plicity, it is convenient to treat a two-dimensional circular condensate that is uni- 
form in the z direction over a length Z. The full condensate wave function W{r,z) 
can then be written as yfN/Z\^{r), where \i/{r) is a two dimensional wave func- 
tion with unit normalization / c/^r | i//p = 1 . 

The general two-dimensional GP energy functional in the rotating frame be- 
comes 



E'[\ir,Q]= fd^rxjf 



( \ 



\ one-body oscillator i^' interaction J 



V, (18) 



where p = — r'^V, L^ = z-r xp, and g2D = Ng/Z. The one -body oscillator hamil- 
tonian J^' in the rotating frame is exactly soluble and has the following eigenval- 

29 

ues^ 

£„,„ =/j[ft)^+n(ft)^+X2)+m(ft)^-X2)], (19) 

where n and m are non-negative integers. In the limit Q. — > 03^, these eigenvalues 
are essentially independent of m, which implies a large degeneracy. The other 
quantum number n then becomes the Landau-level index. The lowest Landau level 
with « = is separated from the higher Landau levels by a gap ^ 2h(0^. 

The large radial expansion means a small central density «(0), so that the in- 
teraction energy g2D«(0) eventually becomes small compared to the gap 2h(0±. 
In this limit, it is natural to focus on the lowest Landau level (LLL), with « = 
and non-negative in > 0. The ground-state wave function is a Gaussian y/oo °^ 
exp(— r^/2c/^), where d^ = y/h/Mmj^ is analogous to the magnetic length in the 
original Landau problem of an electron in a uniform magnetic field. The general 
LLL eigenfunctions have a very simple form 

V^o™(r)ocr''V''"^e-'^/24. (20) 

Since rexp(/(/)) is just the polar form of the complex variable ^ = x+ iy, these 
LLL eigenfunctions become i//o„, <== i^'"exp(— r^/2(i^) with m > 0. Apart from 
the ground-state Gaussian, this is just ^'", a non-negative power of the complex 
variable ^. 

Assume that the GP condensate wave function is a finite linear combination of 
these LLL eigenfunctions 



V^(r)LLL = I c,„Yo,n{r) = fiQe-'^^^''^, (21) 
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where /(i^) =Lm>ofmC'" is an awa/yfic/Mncfion of the complex variable i^. Specif- 
ically, f{Q is a complex polynomial and can be factorized as /(i^) = ]"[/ (C ^ 0) 
apart from an overall constant. Note that f{Q vanishes at each of the points {^j}, 
which are the positions of the nodes of of V'lll- In addition, the phase of the wave 
function increases by 271 whenever ^ moves around any of these zeros {^j} in 
the positive sense. As a result, the LLL trial function (1211) has singly quantized 
vortices located at the positions of the zeros {C/}- 

The mean-field LLL regime implies the striking result that the spatial distri- 
bution of the vortices completely determines the spatial variation of the number 
density n{r) — \ V'lll(''') P- The core size is comparable with the intervortex spac- 
ing / = y^h/MQ,^^ which is the same as d^ in the limit 12 w ft)^. Unlike the 
mean-field TF regime at lower £2, here the wave function i/zLll includes all the ki- 
netic energy. Since the LLL wave functions play a crucial role in the quantum Hall 
effect (two-dimensional electrons in a strong magnetic field), this LLL regime is 
sometimes called the "mean-field quantum Hall" regime. 

It is important to emphasize that we are still in a regime governed by the GP 
equation, so there is still a BEC with a macroscopic condensate wave function. 
The corresponding many-body ground state is simply a Hartree product with each 
particle in the same one-body solution XffhLhi'r): 

N 
fcP (t-1 , r2 , • • • , tat) °= fl l/LLL (Vj). (22) 

y=i 

This is a coherent superfluid state since the GP single -particle state i//lll has 
macroscopic occupation. 



8 Beyond the GP picture: quantum phase transition to highly correlated 
states 

As Q increases still closer toward (0^, the question of what happens beyond the 
mean-field LLL regime remains a subject of vigorous debate. ^^'^^ Generally, a 
quantum phase transition is predicted to take place from the coherent many-body 
BEC ground state in Eq. ( 122b to one of various correlated many-body states that 
are not superfluid and do not have macroscopic occupation. 

To quantify the discussion, it is conventional to define the ratio V = N/N^, 
which is the number of atoms per vortex. Because of similarities to the two- 
dimensional electron gas in a strong magnetic field, the ratio V is called the "filling 
fraction." Current experiments have A^ ^ 10^ and A',, ^ a few hundred vortices so 
that the typical V '^ a few hundred. 

Numerical studies for small number of vortices (A^,, < 8) and variable number 
of particles A' indicate that the GP coherent state is favored for V ^6 — 8.^-^ For 
smaller V, however, the ground state typically has a very different form. Specifi- 
cally, analytical solutions for small A'^ and large angular momentum L are not of 
the factorized GP form. ^^ 

For V < 6 — 8, the ground state is predicted to be one of a sequence of highly 
correlated states similar to some of those known for the quantum Hall effect for 
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electrons. One particularly simple example is a bosonic version of the Laughlin 
state 



Zn\^ 



n<n' \ n=l ^"l 



(23) 



where z„ = x„ + i_y„ refers to the nth particle. The original Laughlin state for elec- 
trons had a power 3 in the double product to ensure antisymmetry, whereas the 
present power 2 ensures symmetry, as appropriate for bosons. These correlated 
many-body states are qualitatively different from the coherent GP form. Specifi- 
cally, the double product in ( 123b involves A^(A'^ —\)/2 factors for all possible pairs 
and vanishes whenever any two particles are close together. This last factor is the 
source of correlations, for it reduces the energy in the typical case of short-range 
repulsive potentials. 

How might one reach the correlated regime? The essential step is to reduce 
the ratio v = N /N^ (the number of atoms per vortex). One possibility is to use an 
elongated condensate with a relatively large vortex array. Subsequent application 
of an optical lattice along the rotation axis would leave an array of thin vortex- 
filled disks that might achieve this goal. 

Another possible idea is to use laser-induced synthetic vector potentials that 
can mimic the effect of rotation. 35,36,37,38 Spielman's group at NIST (Washington 
DC) has indeed produced vortices with this scheme, ^^ although they do not see 
regular arrays. This intriguing approach will certainly receive more attention. 
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